Kramers- Wannier symmetry and strong-weak-coupling duality in the 

two-dimensional ^'^ field model 



B.N.Shalaev 

A.F.Ioffe Physical& Technical Institute, Russian Academy of Sciences, 
Polytechnicheskaya str. 26, 194021 St. Petersburg (Russia), 
Istituto Nazionale di Fisica Nucleare, Struttura di Pavia, Via Bassi, 6-27100, Pavia 

(Italy) 

and Max-Planck Intsitute, Heisenbergstrasse,l, D-70529, Stuttgart, (Germany) 

Abstract 

It is found that the exact beta-function P{g) of the continuous 2D g^^ model 
possesses two types of dual symmetries, these being the Kramers- Wannier (KW) 
duality symmetry and the weak-strong-coupling symmetry f{g), or S-duality. All 
these transformations are explicitly constructed. The /S-duality transformation f(g) 
is shown to connect domains of weak and strong couplings, i.e. above and below g* 
with g* being a fixed point. Basically it means that there is a tempting possibility 
to compute multiloop Feynman diagrams for the /^-function using high-temperature 
lattice expansions. The regular scheme developed is found to be strongly unstable. 
Approximate values of the renormalized coupling constant g* found from duality 
symmetry equations are in good agreement with available numerical results. 
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1 INTRODUCTION 



The 2D Ising model and some other lattice spin models are known to possess the remark- 
able Kramers- Wannier(KW) duality symmetry, playing an important role in statistical 
mechanics, quantum field theory [|l], ^ ^ as well as in the superstring theory 0]. The 
self-duality of the isotropic 2D Ising model means that there exists an exact mapping 
between the high-T and low-T expansions of the partition function 0]. In the transfer- 
matrix language this implies that the transfer-matrix of the model under discussion is 
covariant under the duality transformation. If we assume that the critical point is unique, 
the KW self-duality would yield the exact Curie temperature of the model. This holds 
for a large set of lattice spin models including systems with quenched disorder (for a 
review see [§, |^ ) . Recently, the Kramers- Wannier duality symmetry was extended to the 
continuous 2D g^'^ model ^ in the strong-coupling regime, i.e. for g > g^. 

This beta-function f3{g) is to date known only in the five-loop approximation within 
the framework of conventional perturbation theory at fixed dimension d = 2 ^. 

The strong coupling expansion for the calculation of the beta-function of the 2D scalar 
g^^ theory as an alternative approach to the standard perturbation theory was recently 
developed in @, 0] . 

It is well known from quantum field theory and statistical mechanics that any strong 
coupling expansions are closely connected with the high-temperature (HT) series expan- 
sions for lattice models. From the field-theoretical point of view the HT series are nothing 
but strong coupling expansions for field models, lattices to be considered as a technical 
device to define cut-off field theories (see f^, ^ and references therein) . 

Calculations of beta-functions are of great interest in statistical mechanics and quan- 
tum field theory. The beta-function contains the essential information on the renormalized 
coupling constant g*, this being important for constructing the equation of state of the 
2D Ising model. Duality is known to impose some important constraints on the exact 
beta- function [p]] . 

In this paper we study other duality symmetries of the beta- function l3{g) for the 2D 
(7$^ theory regarded as a non-integrable continuum limit of the exactly solvable 2D Ising 
model. The main purpose is to construct exlicitly the weak-strong (WS) coupling duality 
transformation f{g) connecting domains of weak and strong couplings, i.e. above and 
below g^:. The last transformation allows one to compute unknown yet multiloop orders 
(6,7,...) of the /^-functions on the basis of lattice expansions ([^ ^). 

The paper is organized as follows. In Sect. II we set up basic notations and define 
both the correlation length and beta-function (3{g). In Sect. Ill the duality symmetry 
transformation g = d{g) is derived. Then it is proved that (3{d{g)) = d'{g)(3{g). An 
approximate expression for d{g) is also found. Sect. IV contains an explicit derivation 
of the weak-strong coupling transformation whilst in Sect.V in order to illustrate our 
approach the sixth-order term of (3{g) is approximately computed. Sec. VI. contains 
disussion and some concluding remarks. 
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2 CORRELATION LENGTH AND COUPLING CONSTANT 



We begin by considering the classical Hamiltonian of the 2D Ising model (in the absence of 
an external magnetic field), defined on a square lattice with periodic boundary conditions; 
as usual: 

H = -J (2-1) 

<«,i> 

where < i,j > indicates that the summation is over all nearest-neighboring sites; cxj = ±1 
are spin variables and J is a spin coupling. The standard definition of the spin-pair 
correlation function reads: 



G{R) =< (TRao > (2.2) 

where < ... > stands for a thermal average. 

The statistical mechanics definition of the correlation length is given by [|l^ 

2 dluGjp) 

e - ^^lp=o (2-3) 

The quantity is known to be conveniently expressed in terms of the spherical moments 
of the spin correlation function itself, namely 

/iz = ^(/?/a)'G(R) (2.4) 

R 

with a being some lattice spacing. It is easy to see that 

e = -^ (2.5) 

where d is the spatial dimension (in our case d = 2). 

In order to extend the KW duality symmetry to the continuous field theory we have 
need for a "lattice" model definition of the coupling constant g, equivalent to the con- 
ventional one exploited in the RG approach. The renormalization coupling constant g of 
the g(^^ theory is closely related to the fourth derivative of the "Helmholtz free energy", 
namely d^F{T, ■m)/dm^, with respect to the order parameter m = ($). It may be defined 
as follows (see [Q) 

9iT,h) = -^-^ + 3—-^^ (2.6) 

where x is the homogeneous magnetic susceptibility 

X = J d'xGix) (2.7) 

It is in fact easy to show that g{T, h) in Eq.( p^.6D is merely the standard four-spin corre- 
lation function taken at zero external momenta. The renormalized coupling constant gc 
of the critical theory is defined by the double limit 
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g, = hm\imgiT,h) (2.8) 

ft— *0 T^Tc 

and it is well known that these limits do not commute with each other. As a result, Qc is 
a path-dependent quantity in the thermodynamic (T, h) plane . 



Here we are mainly concerned with the coupling constant on the isochore line g{T > 
Tc, /i = 0) in the disordered phase and with its critical value 

g*= hm o(T,/. = 0) = -^^^|^U=o (2.9) 
Here for convenience we have used a simpler notation g* in contrast to the well-known 



paper where authors made use another symbol g*^. 
3 KRAMERS- WANNIER SYMMETRY 

The standard KW duality tranformation is known to be as follows [jl|, ||, || 

sinh(2i^') = — (3.1) 
^ ^ sinh(2is:) ^ ^ 

We shall see that it will be more convenient to deal with a new variable s = exp{2K) tanh(K), 
where K = J/T. 

It follows from the definition that s transforms as s = 1/s; this implies that the 
correlation length of the 2D Ising model given by = -(jz^ is a self-dual quantity 
Now, on the one hand, we have the formal relation 

= (3.2) 

dg 

where s{g) is defined as the inverse function of g{s), i.e. g{s{g)) = g and the beta-function 
is given, as usual, by 

e§ = m (3.3) 



On the other hand, it is known from 



d^ ^ 2.(1-.) 
^d^ (l + s) ^ ^ 



From Eq.s ( |3.2| ) - ( p.4|) , a useful representation of the beta- function in terms of the 



s{g) function thus follows 

_ 2s{g){l - s{g)) 
^^'^ " (1 + s{g)) {ds{g)/dg) ^^-^^ 

From this representation it follows that l3{s{g)s(g)=i = 1 It is reasonable to assume (but it 
is not obvious!) that this fixed point is a simple zero. Then from this equation it would 
follow that uj = P'{g*) = 1 in agreement with the classical paper |]T3|. 
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Nevertheless, it should be mentioned that this remarkable result " contradicts" to con- 
formal field theory results predicting that for the two-dimensional Ising model = |. In 
the prominent book 0, being the Bible for experts in quantum field theory and critical 
thenomena, it stands just a question mark instead of uj. (Another approach one may find 
in |16[) and discussion in 0). 

Let us define the dual coupling constant g and the duality transformation function 
d{g) as 



s{9) = -Y^; 9 = d{g) = s (3.6) 

where s~^{x) stands for the inverse function of a; = s{g). It is easy to check that a 
further application of the duality map d{g) gives back the original coupling constant, i.e. 
d{d{g)) = g, as it should be. Notice also that the definition of the duality transformation 
given by Eq. (|3.6| ) has a form similiar to the standard KW duality equation, Eq. ( |3.1|) . 

Consider now the symmetry properties of l3{g). We shall see that the KW duality 
symmetry property, Eq. (|3.1|) , results in the beta-function being covariant under the 
operation g d{g): 

f3{d{9)) = d\g)m (3.7) 
To prove it let us evaluate [3{d{g)). Then Eq. (pl5|) yields 

= , '^'f , (3.8) 

(l + s(g))(da(g)ldg) 

Bearing in mind Eq. ( |3.6| ) one is led to 

(^id{g)) = , ^-^/^-^ (3.9) 

s{g){\ + s{g)){ds{g)ldg) 

The derivative in the r.h.s. of Eq. (|3.9| ) should be rewritten in terms of s((y') and d(^g). It 
may be easily done by applying Eq. (p.6|) : 



ds{g) _ d \ _ s'{g) 1 
dg dgs{g) s^{g)d'{g) 

Substituting the r.h.s. of Eq. ( |3.10| ) into Eq. ( |3.9| ) one obtains the desired symmetry 
relation, Eq. (|3.7| ). 

Therefore, the self-duality of the model allows us to determine the fixed point value 
in another way, namely from the duality equation d{g*) = g*. 
Making use of a rough approximation for s{g), one gets [^j 



, , 2 24 2 1 2 , , 

so ~- + — ~ - = 3.11 

9 9^ 9I-IV9 9-12 ^ ^ 

Combining this Pade-approximant with the definition of d{g), Eq. (|3.6| ), one is led to 
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«=4|3# (3-12) 

The fixed point of tliis function, d{g*) = g*, is easily seen to be g* = 14. Tlie recent 
numerical and analytical estimates yield g* = 14.69 (see [13], |T^ and references 

therein) . 

It is worth mentioning that the above-described approach may be regarded as another 
method for evaluating g*, fully equivalent to the standard beta-function method. 

4 STRONG- WEAK-COUPLING DUALITY 

The beta-function of the model under discussion possesses a specific algebraic property 
Eq. (|3.5|) (KW duality) which allows to construct the weak-strong-duality transformation 
f{g) connecting both the weak-coupling and strong coupling regimes. 

Nowadays both the five- loop approximation results (|l3]) and the strong coupling 
expansion for the beta-function |^ are known rather well. These are given by 



p^(g) = 2^ -2/ + 1.432347241/ - 1.861532885c/^ 

+ 3.164776688/ - 6.520837458c/^ + O(/) (4.1) 



, 12 9 27 81 3645 15309 2187 7, ^ 

/32ig) = -2(7+ + + - - + + Oig-') (4.2) 

TT Tv'^g TV'^g'^ QTv^g-^ Ibn^g^ 62'n-°g° bA7r'g° 

Here indices 1,2 stand for the weak and strong coupling regimes respectively. The 
main goal of this Section is to determine a dual transformation f(g) such as f[f{g)] = g 
relating beta-functions Pi{g) and P2{g)- 

From Eq. ( p.5D one can easily find the functions Si{g), S2{g) and their inverse functions 
Gi{s) = S{^{g),G2{s) = S2^{g) corresponding to two regimes. Simple but cumbersome 
calculations lead at 



Gi{s) = s + + 0.35808681045^- 0.1166327797/ 

- 0.1968226859s^ - 0.1299831557s^ + O(s^) 

S^{g) = ^-/ + 1.6419131896/- 3.09293317/ 

+ 6.361881481^^ - 13.78545095^^ + 0{g^) 

s G[0,1]; 9e[0,g*] (4.3) 
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^ , , 3 9 9s 18s2 108s3 

^2 5 = -— + — - — + 



4:713 2n 4:71 n tt 
618s^ 34745^ 194945^ 



+ + + 0{s' 

TT TT TT 

^ , , 2 x 3 24 X 32 264 X 33 2976 x 3^ 

32(9) = -E + ,0 X9 + /o X. + 



4 



^71(7 (871(7)2 (871(7)3 (871(7) 

35136 X 3^ 423680 x 3*^ 5149824 x 3^ 63275520 * 3^ 



(87r(7)5 (87r(7)6 {87,gy {87rgY 

+ 0{g-') 

s G[0,1] 9e[g*,oo) (4.4) 

Having been equipped with these formulas one may easily construct two branches of 
the same duality transformation function fuig) and f2i{g) defined in different domains 
of g. The functions are 

1 1 47ro 287ro2 

- ^ ^ +220.5059303^^ 



f2i{9) G2{S,{g)) 3 3 

- 1766.8145^^ + 14816.94007^^ - 127842.5955/ 
9 e[0,(7*]; /2i(^7)G [(7*,oo] (4.5) 



, , , ^ / o / ^^ 3 63 0.61714739472 

fuig) = Gi{S2{g)) = — + -—r^ + 



4ng ' IQn^g^ g^ 

0.9560453953 1.502156783 2.368311503 
+ 2 + ^ + S 

g^ 

+ 0{g') 

g G[(7*,oo); fM e[Q,g*] (4.6) 

Functions found above look like inversion, but they are not so trivial. An interesting 
nontrivial example of the 2D model disordered Dirac fermions was discovered in [0|. It 



was shown that the beta-function of the (nonintegrable) model under consideration also 
exhibits the strong- weak coupling duality such as g* ^ ^ [ll7| . 

It is worth noting that the transformation found is a self-dual indeed 

/l2(/2l(^7)) = /2l(/l2(^7))=^? (4.7) 

Moreover, by definition weak-strong coupling expansions of j3{g) are related to each 
other in the following way: 

= ^-^^-^ (4.8) 

ft(j) = '^2^ (4.9) 
These equations are basic results of these article. 
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5 HIGHER-ORDER TERMS OF THE BETA-FUNCTION 



It rather amusing that Eq. (|4.6| ) looks hke a geometric series. Making use of the Pade 
method we arrive at 



fi2ig) 



0.2387324146^2 _ 0.0745907136^ + 0.0850867165 
- 1. 983571 753C/2 + 1.086109562^ - 0.6919672492 

9 e[g*,ooy, fi2{9)e[0,g*] (5.1) 



The weak-strong-couphng duahty equation and strong-couphng expansion yield the 
following numerical values 



fi2{g*)-9* = 0; g* = U.38 

P2{g*) = 0; ^7* = 14.63 (5.2) 

being in good agreement with modern estimates [|18], 0, 

Finally, let's consider how one can compute the l3{g) in the multiloop approxima- 
tion via the strong-coupling expansion and the S-duality function. In order to find that 
one should exploit Eq.([1.8|),Eq.(^^) and the approximate expression for fi2{g) given by 
Eq.(U. 

After some tedious but routine calculations we arrive to some polynomial of 7th degree 
for /3,{g). 



(3^(g) = 2^ - 2/ + 1.432347241c/^ - 1.861532885^^ 

+ 3.164776688/ - 6.520837458/ - 331.454743/ (5.3) 

It is easily seen that the first 6 terms excepting for the 7th one are the exact per- 
turbation expansion for Pi{g) [|10|- It would be tempting but wrong to regard Eq. (|573|) 
as a /3((?)-function in the 7th loop approximation. In fact, the function in Eq. ( |5.1| ) is 
approximate, so that we have to estimate an accuracy of our calculations. 

Let's illustrate our treatment by using a simple example. Suppose, that a difference 
between the "exact" duality function /ff"'^*(5') and the approximate one given by Eq. (|5?T|) 
reads 



^^act( 0.2387324146/ - 0.0745907136^ + 0.0850867165 ^ 

•^^2 ~ - 1.983571753/ + 1.086109562^ - 0.6919672492 ^ ^ ' ^ 



with b being an arbitrary unknown parameter. The straightforward calculation shows 
that a "new" 7th loop contribution computed by making use of the Eq. (|5.4| ) depends on 
the fitting parameter b and differs vastly from the previous one, it being 



f3^(g) = 2^- 2/ + 1.432347241/- 1.861532885/ 

+ 3.164776688^^ - 6.520837458/ + (-331.454743 + 271519.8038076)^^(5.5) 
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Assuming for instance that b = 10~*, so that the difference between exact dual function 
/i2"'^*(fl') ^^'i ^ thai one is numerically small indeed, one obtains that a change of the sixth 
loop coefficient is of order 1. 

Thus, we see that the approach suggested above provides a regular scheme for com- 
puting higher-order corrections to the /5((7)-function on the basis of lattice high-order 
expansions. In other words, one obtaines a tempting possibility to compute (approxi- 
mately) multiloop Feynman diagrams on the basis of Eg. ( [4 .91 ) and of high-temperature 
expansions [Q. The serious drawback of that scheme is that it being unstable from a 
mathematical point of view. 

6 CONCLUDING REMARKS 

We have shown that the /3-function of the 2D g^^ theory does have two types of dual 
symmetries, these being the Kramers- Wannier symmetry and the weak-strong-coupling 
symmetry or S-duality. These important symmetries were shown to have a thorough 
mathematical background, it being highly nontrivial algebraic structure Eg. ( ^.51 ) of l3{g). 

Our proof of the KW symmetry is based on the properties of g{s), s{g) defined only for 
l<s<oo;(7*<5f<oo and therefore doesn't cover the weak-coupling region, < g < g*. 
So that the statement is that the beta-function P{g) possesses the KW symmetry only in 
the strong-coupling region. 

In contrast to widely held views, the KW symmetry imposes only mild restrictions 
on f3{g). It means that this symmetry property fixes only even derivatives of the beta- 
function (3'^'^^\g*)(k = 0,1,...) at the fixed point, leaving the odd derivatives free, in 
particular, the critical exponent u, responsible for corrections to scaling. 

We have established the existence of the nontrivial weak-strong-coupling dual function 
f{g) connecting two domains of both weak coupling and strong coupling. Given both 
perturbative RG calculations and lattice high-temperature expansions that S-function 
f{g) can be approximately computed. We also explicitly computed high-order terms of 
P{g). In summary, the value of the above developed approach is considerably deduced by 
its being strongly unstable. 
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